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straight  duct  was  tested  with  an  acoustic  liner  consisting  of  a  matrix  of 
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the  hard  walled  configurations  while  there  were  some  discrepancies  with  the 
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wall  results  in  mainly  due  to  the  particular  liner  theory  used  to  calculate 
the  effective  admittance  of  the  liner. 


UNCLASSIFIED 


secuaiTv  ctMtincATio*  o* 


Abstract 


This  report  summarizes  the  work  performed  in  the  third  and  last  year 
of  an  AFOSR  sponsored  research  program  (AFOSR  Contract  No.  F49620-77- 
C-0066).  This  research  program  was  concerned  with  the  development  of  an 
analytical  technique,  based  on  an  integral  representation  of  the  external 
solutions  of  the  Helmholtz  equation,  for  the  prediction  of  the  sound  radiated 
from  complicated,  acoustically  lined,  axisymmetric  bodies  having  complex 
sound  sources.  The  purpose  of  this  research  program  was  to  generate 
efficient  computer  codes  for  the  prediction  of  the  sound  radiated  from 
acoustically  lined  jet  engine  inlets.  During  the  first  two  years  of  work  under 
this  contract  the  above  goals  were  accomplished  and  are  documented  in  the 
previous  two  AFOSR  ANNUAL  TECHNICAL  REPORTS  numbered  AFOSR- 
TR-7S-0696  and  AFOSR-TR-79-0614.^  _f\  Q  S  ‘  > 

This  report  is  concerned  with  the  progress  made  during  the  third  year 
of  this  contract  when  experimental  tests  were  run  on  various  configurations 
for  comparison  with  and  verification  of  the  results  obtained  from  the 
axisymmetric  computer  codes  which  model  the  integral  technique.  In  the 
experimental  tests  two  geometrical  configurations  were  studied,  a  straight 
duct  and  a  jet  engine  inlet.  Both  of  these  configurations  were  tested  with 
hard  walls  and  the  straight  duct  was  tested  with  an  acoustic  liner  consisting 
of  a  matrix  of  Helmholtz  resonators.  It  was  found  that  very  good  agreement 
was  obtained  for  the  hard  walled  configurations  while  there  were  some 
discrepancies  with  the  lined  wall  case.  It  is  conjectured  that  this  discrepancy 
in  some  of  the  lined  wall  results  is  mainly  due  to  the  particular  liner  theory 
used  to  calculate  the  effective  admittance  of  the  liner. 


Introduction 

This  report  summarizes  the  resuits  obtained  during  the  third  year  of 
support  under  AFOSR  contract  number  F49620-77-C-0066.  This  contract  was 
initiated  February  1,  1977  and  the  results  obtained  during  the  first  two  years 
of  support  are  contained  in  AFOSR  technical  reports  AFOSR-TR-78-0696 
and  AFOSR-TR-79-0614. 

The  main  objective  of  the  research  program  conducted  under  this 

contract  was  to  develop  an  analytical  technique,  both  the  theory  and 

associated  computer  codes,  for  predicting  the  sound  field  radiated  from 

axisymmetric  jet  engine  inlet  configurations  with  lined  walls  and  to  compare 

some  analytical  predictions  with  the  results  of  experimental  tests.  The 

development  of  the  theory,  which  is  based  on  a  special  integral 

representation  of  the  external  solutions  of  the  Helmholtz  equation,  was 

motivated  by  the  need  for  an  analytical  approach  that  could  be  used  to 

predict  the  effects  of  sound  source  modifications  and  acoustic  liners  on  the 

sound  field  radiated  from  an  inlet  without  having  to  resort  to  costly,  full 

scale  experimental  testing.  During  the  first  two  years  of  this  contract  the 

(1  3) 

axisymmetric  formulation  of  the  integral  solution  technique  ’  was 
developed  along  with  two  efficient,  general  computer  programs;  one  to  solve 
for  the  surface  distributions  of  the  acoustic  quantities  of  interest  and  the 
other  to  solve  for  the  distributions  of  the  acoustic  quantities  in  the  field 
surrounding  the  body. 

During  the  third  contract  year  sound  radiation  experiments  were 
performed  with  two  geometrically  different  configurations,  a  straight  duct 
and  the  QCSEE  jet  engine  inlet  of  Ref.  (4),  for  comparison  with  analytical 
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results  for  the  same  bodies.  The  straight  duct  configuration  was  tested  with 
both  a  hard  and  a  soft  wall  while  the  inlet  configuration  was  only  tested  with 
a  hard  wall.  Tasks  performed  during  the  third  contract  year  included: 

A. )  Calibration  of  the  anechoic  chamber 

This  was  necessary  to  find  out  how  well  the  anechoic  chamber 
approximated  the  results  that  would  be  gotten  if  true  free  field 
measurements  could  be  made.  These  calibrations  were  done  with  all  of  the 
support  apparatus  for  the  test  bodies  and  the  microphones  in  place  so  as  to 
get  an  idea  of  the  magnitude  of  the  errors  involved. 

B. )  Set-up  of  the  electronic  equipment 

This  task  was  not  as  straight  forward  as  it  might  seem  as  it  was  found 
that  in  order  to  get  relatively  stable  amplitude  and  phase  measurements  it 
was  necessary  to  pre-condition  the  signals  coming  from  the  microphones 
through  the  use  of  both  high  and  low  pass  filters.  Since  all  of  the  signals 
were  passed  through  the  same  filters,  including  the  reference  signal  from 
the  reference  microphone,  any  phase  or  amplitude  shift  caused  by  the  filters 
could  be  subtracted  out. 

C. )  Conduct  of  the  experiments 

Each  of  the  experimental  tests  were  run  at  least  twice  on  two 
different  days  and  the  results  of  these  separate  tests  were  compared.  This 
was  necessary  to  get  an  estimate  of  the  repeatability  of  the  experimental 
data.  It  was  also  used  as  a  check  on  the  validity  of  the  data  as  the  data  was 
read  directly  from  meters  and  transcribed  by  hand  as  no  on  line  data 
acquisition  system  was  available  for  use  on  this  project. 
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D. )  Prediction  of  the  tested  liner  admittance 

A  liner  was  designed  and  fabricated  for  the  lined  straight  duct  tests. 
To  do  this  design  the  theory  of  Garrison^  was  employed.  The  liner  thus 
designed  was  a  matrix  of  Helmholtz  resonators  which  had  its  absorption  peak 
below  the  IT  mode  of  the  duct  as  this  was  where  all  of  the  testing  was  done. 

E. )  Data  Reduction  and  comparison  with  theoretical  predictions 

The  experimental  data  was  both  taken  and  reduced  by  hand.  The  data 
reduction  was  not  a  very  time  consuming  process  as  all  that  had  to  be  done 
was  subtract  the  reference  microphones  values  of  the  amplitude  and  phase 
from  the  values  measured  at  the  driver  plane  and  the  radiation 
measurements  in  the  field.  This  was  done  so  as  to  negate  any  shift  in  the 
measured  amplitude  or  phase  by  the  signal  filters. 

The  computer  programs  were  then  run  for  the  same  conditions  (i.e., 
driver  power  and  frequency)  as  the  experimental  tests.  These  results  were 
then  compared  with  the  experimental  results  and  a  brief  error  analysis  was 
done. 

F. )  Determination  of  the  "effective"  liner  admittance 

The  results  of  the  experimental  tests  and  the  theory  were  found  to  be 
in  very  good  agreement  except  near  the  calculated  absorption  peak  of  the 
liner.  It  was  determined  that  the  Garrison  theory^  used  to  calculate  the 
theoretical  liner  admittances  was  probably  at  fault.  Systematic  computer 
runs  were  then  performed  parametrically  varying  the  admittance  until  good 
agreement  was  obtained  between  the  theoretical  and  experimental  results. 
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Included  in  this  report  is  a  list  of  all  the  publications  generated  by 
this  research  effort.  This  is  included  as  Appendix  A.  Also,  Appendix  B 
contains  a  list  of  all  the  conference  presentations  pertaining  to  the  work 
performed  under  this  contract  and,  finally,  a  copy  of  the  most  recent 
publication  that  was  a  direct  result  of  the  research  performed  in  the  past 
year  on  this  contract  is  included  as  Appendix  C. 
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A.  Calibration  of  the  Anechoic  Chamber 

The  anechoic  chamber  is  10'  by  13'  by  6  1/2'  in  height.  It  has 
approximately  2'  of  sound  absorbing  fiberglass  insulation  in  the  floor,  walls, 
and  ceiling.  A  plan  view  of  the  anechoic  chamber  is  presented  in  Fig.  8  of 
Appendix  C. 

It  was  necessary  to  calibrate  the  anechoic  chamber  with  all  of  the 
microphone  stands  and  the  mounting  stand  for  the  test  apparatus  in  place  so 
that  an  estimate  of  the  experimental  errors  in  the  free  field  measurements 
caused  by  this  apparatus  could  be  made.  A  University  acoustic  driver  was 
mounted  vertically  on  the  support  pedistal  for  the  test  apparatus,  37.5"  off 
the  floor,  so  that  it  would  radiate  sound  equally  in  all  directions,  thus 
approximating  a  simple  source.  A  Vi '  diameter  Briiel  and  Kjaer  microphone 
was  placed  42.5"  from  the  driver  in  the  field  and  was  used  as  a  reference 
microphone  for  the  amplitude  readings. 

A  K"  diameter  B  <5c  K  microphone  was  then  moved  in  the  direction 
that  the  test  ducts  would  face,  to  be  henceforth  known  as  direction  1,  (See 
Fig.  8,  Appendix  C  and  Fig.  1.)  and  measurements  were  taken  every  5"  out  to 
45".  This  procedure  was  repeated  in  the  direction  of  the  door  of  the  chamber 
(90°  to  the  original  direction  of  travel)  and  again  measurements  were  taken 
every  5"  up  to  45"  from  the  driver.  This  will  be  referred  to  as  direction  2.  At 
each  microphone  position  measurements  were  made  of  the  amplitude  and 
phase,  relative  to  the  driver  input  signal,  at  frequencies  from  300  to  750  Hz 
at  50  Hz  increments  as  the  tests  were  all  to  be  run  in  this  range.  It  should  be 
noted  here  that  the  lower  limit  was  imposed  by  the  University  driver  and 
that  the  upper  limit  was  imposed  by  the  requirement  that  all  of  the  tests  be 


Test  Duct  Position  (Not  Present 
During  the  Calibration  Runs) 

_ Driver  Position  for 
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Fig.  1.  Plan  View  of  Anechoic  Chamber  with  Driver  and 
Microphone  Positions  for  the  Calibration  Runs. 
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run  below  the  IT  mode  of  the  inlet  and  the  straight  duct.  The  latter 
requirement  was  to  facilitate  the  data  reduction  and  comparison  with  the 
theoretical  model  and  the  fact  that  we  had  no  simple  method  of  driving 
higher  modes  in  the  test  ducts. 

Results  of  the  amplitude  calibration  in  direction  1  are  presented  in 
Fig.  2.  The  solid  lines  are  the  theoretical  curves  for  the  Sound  Pressure 
Level  in  decibels  for  a  simple  source  in  a  free  field.  The  equation  for  these 
curves  is 

SPL(dB)  =  -  20  Log  r 

(1) 

where  r  is  the  distance  from  the  source.  In  this  plot  the  amplitude  is 
referenced  to  the  5"  position.  The  worst  point  in  this  direction  is  at  30"  and 
350  Hz  where  the  measured  reading  is  3.5  dB  below  the  theoretical  exact 
curve.  All  of  the  test  runs  were  made  with  the  microphones  at  40"  where 
the  error  was  always  less  than  3  dB  at  all  frequencies. 

The  results  of  the  amplitude  calibration  in  direction  2  are  presented 
in  Fig.  3.  The  worst  point  in  this  direction  occurs  at  30"  and  400  Hz  where 
the  measured  sound  pressure  level  is  9.5  dB  below  the  free  field  prediction. 
Again,  all  of  the  test  runs  were  made  with  the  microphones  at  40"  where  the 
error  was  always  less  than  3  dB. 

In  Figs.  4  and  5  the  phase  comparisons  are  presented.  The  X's  in  these 
plots  denote  points  where  no  stable  phase  reading  could  be  obtained.  This 
was  later  found  to  be  caused  by  widely  different  signal  strengths  (i.e., 
voltages)  being  fed  into  the  phase  meter.  This  was  corrected  in  subsequent 
tests  by  using  another  microphone  signal  as  the  reference  input  rather  than 
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the  output  from  the  oscillator  that  drove  the  acoustic  driver.  As  can  be  seen 
most  of  the  results  are  within  10°  of  the  values  that  were  calculated  using 
the  simple  source  in  a  free  field  assumption.  Some  values  are  in  error  by  as 
much  as  20°  (See  Fig.  k,  direction  1  at  3( J'  and  35"  from  the  driver  for  300 
Hz.);  however,  these  "bad"  results  are  suspect  as  the  phase  meter  could  not 
be  stabilized  at  adjacent  points.  The  test  results  in  these  plots  are 
referenced  to  the  5"  position  as  it  was  closest  to  the  driver  and  therefore 
subject  to  the  least  interference.  This  is  a  source  of  error  that  could 
potentially  bias  the  results  at  a  given  frequency  since  if  the  phase  reading  at 
the  5"  position  was  high  the  rest  of  the  values  at  all  other  distances  would  be 
low  by  that  amount  and  vice  versa.  It  is  conjectured  that  this  is  what 
happened  at  a  couple  of  frequencies  such  as  600  and  700  Hz  in  Fig.  5. 

From  the  calibration  runs  that  were  performed  in  the  anechoic 
chamber  with  all  of  the  support  equipment  for  the  model  and  the 
microphones  in  place  it  was  concluded  that  the  results  obtained  from  tests  in 
the  chamber  would  be  accurate  to  at  least  3  dB  in  amplitude  and  10°  in 
phase.  As  a  matter  of  fact,  the  results  obtained  from  tests  in  the  chamber 
were  in  general  found  to  be  in  better  agreement  with  the  theoretically 
calculated  "exact"  results  than  the  calibration  runs  would  lead  one  to 
expect. 
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B.  Set-up  of  the  Electronic  Equipment 


The  equipment  set-up  can  be  basically  broken  down  into  two  separate 
entities;  the  sound  generating  equipment  and  the  sound  measuring 
equipment.  During  the  calibration  of  the  anechoic  chamber  (See  previous 
section.)  these  two  systems  were  linked  in  that  the  oscillator  that  produced 
the  driver  signal  was  also  used  as  a  phase  reference  for  the  sound  readings. 
This  arrangement  was  changed  in  all  subsequent  tests  as  it  was  found  that 
due  to  the  large  voltage  difference  in  the  signal  strengths  between  a 
microphone  and  the  oscillator  the  phase  meter  would  not  stabilize  if  the 
signals  were  close  to  180°  out  of  phase.  This  being  the  case,  another 
microphone  was  used  as  the  reference  in  subsequent  tests. 

In  all  of  the  tests,  a  75  watt  University  Sound  heavy-duty  driver  was 
employed.  The  University  driver  was  powered  by  a  Krohn-Hite,  Model  DCA- 
50R,  wide-band  50  watt  amplifier  which  is  capable  of  amplifying  signals 
from  D.C.  to  500  KHz.  A  Model  MT-56  Krohn-Hite  matching  transformer 
was  used  to  match  the  impedances  of  the  amplifier  and  the  driver  to  get 
maximum  power  out  of  the  driver.  The  amplifier  was  triggered  by  a 
Hewlett-Packard  Model  202C  low  frequency  oscillator  which  can  produce 
signals  from  1  Hz  to  100,000  Hz.  The  frequency  output  of  the  oscillator  was 
constantly  monitored  by  two  time  base  counters;  one  with  fast  response  for 
ease  of  setting  the  desired  frequency  and  one  with  slow  response  for 
checking  the  accuracy  of  the  set  frequency.  The  first  was  a  Hewlett- 
Packard  Model  5302  A  50  MHz  universal  counter  and  the  second  was  a 
Monsanto  Model  104  A  preset/ variable  time-base  counter.  In  all  of  the 
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conducted  tests  the  frequency  was  always  set  to  within  +  0.2  Hz  of  the 
desired  frequency  before  any  measurements  were  made. 

For  the  actual  sound  measurements  two  sizes  of  Briiel  &  Kjaer 
microphones  were  used;  ft  inch  microphones  were  used  outside  the  test  duct 
for  the  field  measurements  and  &  inch  microphones  were  used  inside  the 
duct  for  reference  and  driver  power  measurements  (See  Figs.  4  and  7  of 
Appendix  C.).  Also,  two  different  types  of  pre-amplifiers  were  used  on  these 
microphones;  tube  type  which  require  a  heating  element  voltage  input  and 
transistor  type  which  do  not. 

The  ft  inch  microphones  were  used  with  the  tube  type  pre-amplifiers 
and  the  signals  from  these  microphones  were  fed  into  a  B  &  K  two  channel 
microphone  selector,  Type  4408,  and  from  there  into  a  B  &  K  microphone 
amplifier.  Type  2604.  The  ft  inch  microphones  were  used  with  the  transistor 
type  pre-amplifiers.  In  each  of  the  tests  two  ft  inch  microphones  were  used 
and  five  ft  inch  microphones  were  used.  To  multiplex  the  five  ft  inch 
microphone  signals  an  in  house  fabricated  five  channel  microphone  selector 
was  used  and  these  were  fed  into  a  B  &  K  Type  2606  measuring  amplifier. 
Both  of  the  B  <5c  K  amplifiers  gave  direct  decibel  readings  for  the 
microphones  to  0.1  dB,  although  the  readings  were  certainly  not  reliable  to 
this  accuracy  due  to  amplifier  and  pre-amplifier  drift  (this  will  be  discussed 
in  detail  in  the  following  section).  The  two  microphone  multiplexers  were 
necessary  as  each  of  the  measuring  amplifiers  could  only  handle  one 
microphone  at  a  time. 

Although  the  signals  were  good  enough  to  give  reliable  decibel 
readings  they  had  to  be  pre-conditioned  before  good,  steady  phase 
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measurements  could  be  made.  To  achieve  this,  two  Krohn-Hite  Model  3202  R 
filters  were  used.  Both  of  the  microphone  signals,  one  being  used  as  the 
reference  signal,  were  passed  through  both  a  high  and  a  low  pass  filter. 
These  filters  were  set  for  all  of  the  tests  at  250  Hz  for  the  high  pass  and  800 
Hz  for  the  low  pass  as  all  of  the  tests  were  run  between  300  and  750  Hz.  The 
phase  was  measured  by  a  Wavetek  Model  740  phase  meter  and  was  displayed 
in  degrees  (-180°  to  180°)  on  a  Fluke  8000  A  digital  multimeter. 

Finally,  all  of  the  signals  were  run  through  an  H-P  180  A  dual  beam 
oscilloscope.  This  was  used  to  keep  a  check  on  the  equipment  by  constantly 
checking  the  signal  quality.  Also,  this  was  used  to  check  if  any  of  the 
equipment  was  being  over  driven  or  saturated  at  any  given  set  of  test 
conditions. 


C.  Conduct  of  the  Experiments 

As  stated  before,  each  of  the  tests  was  performed  twice  on  separate 
days  for  verification  of  the  data  as  it  was  read  directly  from  meters  and 
transcribed  by  hand.  The  results  that  are  shown  for  the  experimental  tests 
are  in  fact  an  average  of  the  two  tests  run  for  each  experimental 
configuration.  Also,  the  ambient  temperature  was  monitored  during  each 
test  so  that  corrections  could  be  made  for  the  varying  speed  of  sound  and 
the  ambient  pressure  was  measured  before  and  after  each  run  so  that  the 
characteristic  impedance  of  the  medium  (i.e.,  the  air  in  the  anechoic 
chamber)  could  be  calculated.  The  latter  measurement  was  only  important 
when  a  lined  configuration  was  being  tested. 

Since  nine  measurements  were  made  in  the  field  for  each  test 
frequency  and  since  only  five  Vt  inch  microphones  and  stands  were  available, 
four  of  the  microphones  had  to  be  moved  during  each  test.  The  fifth 
microphone,  on  the  centerline  of  the  test  duct,  was  never  moved  and  was 
used  as  a  reference  to  check  that  the  test  conditions  (i.e.,  the  driver  power 
output)  were  the  same  with  the  other  microphones  in  both  positions.  A 
second  reference  microphone,  a  14  inch  one,  was  placed  in  the  nozzle  section 
of  the  test  set-up  (See  Figs.  4  and  7  of  Appendix  C.)  just  in  front  of  the 
University  driver  so  that  a  second  independent  check  of  the  driver  power 
output  could  be  made  for  each  test  condition. 

The  sequence  of  events  for  each  test  was  the  same  so,  for  the  sake  of 
simplicity,  the  straight  duct  will  be  used  as  the  example  here.  The  only 
difference  lies  in  the  fact  that  the  straight  duct  was  run  with  both  hard  and 
soft  walls  while  the  inlet  configuration  was  only  run  with  a  hard  wall. 
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First,  the  test  duct  was  set-up  in  its  soft  wailed  configuration  (i.e., 
with  its  liner  exposed).  The  microphones  were  then  set-up  at  40  inches  from 
the  center  of  the  entrance  plane  of  the  duct  at  22.5°  increments  from  the 
centerline  of  the  duct  (See  Fig.  8  Appendix  C.).  This  was  accomplished  by 
triangulation  and  the  Yi  inch  microphones  were  placed  within  at  least  one 
half  inch  of  their  ideal  location.  The  microphones  were  then  calibrated  using 
a  B  &  K  pistonphone  acoustic  calibrator  and  the  barometric  pressure  was 
noted.  Then,  the  temperature  in  the  chamber  was  read  and  the  chamber  was 
closed.  Next,  the  oscillator  was  set  to  the  desired  frequency  and  the  driver 
power  was  set  by  making  the  microphone  at  the  imaginary  driver  plane  read 
some  pre-selected  dB  level.  All  of  the  other  microphones  phases  and 
amplitudes  were  then  read  and  recorded.  The  chamber  was  then  opened  and 
the  temperature  was  rechecked.  The  chamber  was  closed  and  the  test  was 
repeated  at  a  new  frequency.  The  inlet  was  run  at  frequencies  from  300  to 
750  Hz  at  50  Hz  increments  while  the  straight  duct  was  tested  over  a 
frequency  range  up  to  700  Hz.  This  was  due  to  the  different  characteristic 
lengths  of  the  two  test  bodies. 

Once  the  desired  frequency  range  had  been  tested,  the  liner  in  the 
duct  was  taped  over  so  that  it  was  no  longer  exposed,  the  microphones  were 
recalibrated,  and  the  whole  test  procedure  was  repeated.  Having  done  this 
the  four  microphones  in  the  field  were  moved  to  their  intermediate 
positions,  recalibrated,  and  another  run  was  made  of  all  the  frequencies. 
Finally,  the  tape  was  removed  from  the  liner,  the  microphones  recalibrated 
again  and  the  final  set  of  measurements  were  made.  Having  done  ail  this, 
the  microphones  were  calibrated  one  last  time  and  the  barometric  pressure 
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was  re-read. 


Each  one  of  these  tests  took  an  entire  day  and  required  only  one 
person  most  of  the  time.  Two  people  were  required  when  the  microphones 
were  calibrated,  one  in  the  chamber  with  the  pistonphone  and  the  other 
outside  recording  the  dB  levels. 
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D.  Prediction  of  the  Tested  Liner  Admittance 

For  the  lined  walled  tests,  a  liner  was  designed  using  the  theory  of 
Garrison.^  The  liner  consisted  of  a  matrix  of  180  Helmholtz  resonators  (See 
Fig.  5  Appendix  C.),  9  axial  rows  by  20  radial  rows.  It  was  designed  so  that 
its  maximum  effectiveness  (i.e.,  resonance  peak)  occurred  within  the  range 
of  frequencies  used  in  testing;  that  is,  between  300  and  700  Hz,  below  the  IT 
mode  of  the  duct  and  above  the  driving  floor  of  the  University  driver.  The 
tests  were  also  all  run  in  the  liners  linear  regime  as  calculated  by  the 
theory/^ 

In  the  Garrison  theory  the  linear  regieme  is  defined  as  where  the 
maximum  acoustic  velocity  at  the  orifice  of  the  Helmholtz  resonator  is 
below  60  ft./sec..  Using  this  criterion  it  is  found  that  the  linear  regime  is 
below  the  decibel  levels  in  Table  I  for  a  sea  level  standard  atmosphere  and 
the  respective  driving  frequencies.  Since  in  all  of  the  runs  the  decibel  level 
was  always  less  than  135  dB  the  linear  liner  theory  was  used. 

Using  the  linear  theory,  the  specific  acoustic  impedance  of  the  liner 
is  given  by 

Z  =  0  -  iX 

(2) 

where  the  specific  acoustic  resistance  is  given  by 


0  = 


4 

ape 


(rrTLpf)%  (1  +  t/d) 


(3) 


See  Nomenclature  at  the  end  of  this  section. 
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Table  I 


Maximum  SPL(dB)  where  a  Helmholtz  resonator  exhibits  linear  behavior 


Frequency  SPL 

(Hz)  (dB) 

300  163 

350  160 

400  157 

450  153 

500  148 

550  143 

600  148 

650  152 

700  155 
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and  the  specific  acoustic  reactance  by 


x  -  2TTf°Ieff  (  _I_  1° 

x  c  a  '  fa  "  f  ' 


In  these  equations  the  effective  orifice  length  is  given  by 


1  «  =  t  +  0.85  d  (  1  -  0.7  V5~  ) 
elt 


and  the  resonant  frequency  of  a  Helmholtz  resonator  by 


For  the  test  conditions  these  are  found  to  be  given  as 


(5) 


(6) 


e  =  o.oi398  VT 


X  =  10.186  - f-  ) 

O 


(7) 


where  the  resonance  peak  of  the  liner  f  is  found  to  be  at  558  Hz.  These 

o 

values  are  related  to  the  acoustic  admittance  of  Appendix  C  defined  as  y  = 
V/S  by 


y  =  -ik<  \  ) 

Z  (8) 

where  Z  is  defined  with  an  inward  facing  normal  and  y  is  defined  with  an 
outward  facing  normal;  thus  the  minus  sign. 
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Using  the  definition  of  the  absorption  coefficient  a ;  that  is, 

4e 

a  -  - * - - - 

(i+e)2  +  x2  (9) 

it  is  found  that  the  liner  absorption  curve  is  very  "peaky"  (See  Fig.  6.).  This 
being  the  case  normal  machining  errors  of  0.0005"  were  introduced  into  the 
admittance  calculation  and  the  resonance  frequency  was  found  to  change  by 
1.5  Hz.  It  was  found  that  this  change  in  the  resonance  frequency  could 
change  the  results  at  550  Hz  (i.e.  the  test  frequency  closest  to  the  resonant 
frequency)  by  3  dB  and  5  degrees. 

This  will  be  discussed  further  in  the  ensuing  sections  as  good 
agreement  was  obtained  between  the  theory  and  experiment  except  at  550 
Hz  for  the  soft  walled  duct.  It  is  felt  that  this  disagreement  is  more  the 
fault  of  the  liner  design  approach  at  frequencies  close  to  resonance  rather 
than  inaccuracies  in  the  integral  equation  formulation. 
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Fig.  6.  Absorption  Coefficient  vs.  Frequency  for  Liner  at 
Test  Conditions. 


100 


Nomenclature 


a  characteristic  length  of  the  body 

c  speed  of  sound 

d  Helmholtz  resonator  orifice  diameter 

f  frequency 

f  resonant  frequency  of  a  Helmholtz  resonator 

i  v^r 

k  wave  number 

L  backing  depth  of  a  Helmholtz  resonator 

*eff  effective  orifice  length  of  a  Helmholtz  resonator 

t  orifice  length  of  a  Helmholtz  resonator 

V  cylindricaliy  symmetric  normal  acoustic  velocity 

y  acoustic  admittance 

Z  specific  acoustic  impedance  of  the  liner 

a  absorption  coefficient  of  the  liner 

a  open  area  ratio  of  the  liner 

a  open  area  ratio  of  a  Helmholtz  resonator 

p  density  of  the  medium 

jl  coefficient  of  viscosity  of  the  medium 

I  cylindricaliy  symmetric  acoustic  potential 

9  specific  acoustic  resistance 

X  specific  acoustic  reactance 
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E.)  Data  Reduction  and  Comparison  with  Theoretical  Predictions 

In  the  data  reduction  the  speed  of  sound,  the  density,  and  the 
coefficient  of  viscosity  were  corrected  for  atmospheric  conditions.  The 
speed  of  sound  determines  the  wave  number  from  the  frequency  and  is 
therefore  extremely  important  in  the  data  reduction  for  all  cases  while  the 
density  and  coefficient  of  viscosity  are  only  important  in  the  theoretical 
admittance  calculation  (See  Eq.  (3).). 

It  should  be  pointed  out  that  the  theoretical  and  experimental  models 
used  in  this  study  were  not  identical  (See  Figs.  3  and  4  for  the  straight  duct 
configurations  employed  and  Figs.  6  and  7  for  the  inlet  configurations  used 
in  Appendix  C.).  The  theoretical  bodies  were  given  hemispherical  rear 
terminations  in  the  interest  of  conserving  computing  time  and  space.  It  has 
been  found,  however,  through  theoretical  studies  using  the  integral  equation 
formulation  of  the  problem,  that  the  exact  form  of  the  rear  termination  of  a 
body  has  little  effect  on  the  sound  field  radiated  to  the  forward  half  plane. 
Since  the  forward  half  plane  is  where  all  of  the  experimental  measurements 
were  made  this  is  not  considered  to  be  a  source  of  major  errors. 

In  both  experimental  models  the  University  driver  was  placed  at  the 
throat  of  a  nozzle  section  (See  Figs.  4  and  7  Appendix  C.)  so  that  the  sound 
waves  are  plane  at  the  driver  plane.  This  was  checked  experimentally  and 
the  sound  waves  at  the  driver  plane  were  found  to  be  plane  within  1  dB  in 
amplitude  and  5°  in  phase.  The  amplitude  was  measured  at  this  plane  and 
was  used  as  the  input  for  the  computer  programs.  The  phase  was  not 
measured  directly  here  as  all  of  the  other  phase  measurements  were 
referenced  to  this  microphone. 
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The  straight  duct  was  tested  in  both  a  hard  walled  configuration  and 
in  a  lined  configuration.  For  the  hard  wailed  tests  the  orifices  of  the 
Helmholtz  resonators  of  the  liner  were  simply  covered  with  tape.  The  inlet 
tested  had  an  L/a  (length/radius)  ratio  of  2.0.  The  exact  mathematical  forms 
of  the  curves  that  make  up  the  inlet  contours  can  be  found  in  Ref.  (4).  It 
should  be  noted  here  that  in  the  experimental  inlet  model  the  centerbody 
was  held  in  place  by  four  small  wing  sectional  struts  set  at  90°  from  each 
other.  Measurements  of  the  driver  amplitude  were  taken  between  two  of 
these  struts. 

A  comparison  between  the  predicted  and  measured  amplitudes  for  the 
hard  walled  inlet  are  presented  in  tabular  form  in  Table  II.  It  should  be 
recalled  that  the  presented  experimental  data  is  actually  the  average  of 
results  obtained  in  two  completely  separate  tests.  The  average  absolute 
errors  between  the  experimental  and  theoretical  results  are  included  at  each 
frequency.  The  results  are  presented  on  a  quarter  circle  with  a  radius  of  40" 
centered  at  the  duct  entrance  plane  (See  Fig.  8  Appendix  C.).  The 
experimentally  measured  amplitude  at  the  driver  plane  is  also  included.  The 
phase  results  are  presented  in  Table  III.  The  driver  phase  is  taken  to  be  zero 
as  it  is  the  reference  phase.  The  results  for  the  inlet  are  plotted  for  600  Hz 
in  Fig.  11  of  Appendix  C. 

As  can  be  seen  from  Tables  II  and  III  the  average  absolute  error 
between  the  theoretical  and  experimental  results  for  the  inlet  configuration 
is  always  less  than  3  dB  in  amplitude  and  10°  in  phase.  At  this  juncture  it 
seems  appropriate  to  discuss  some  of  the  sources  of  error  and  their 
estimated  magnitude  so  as  to  put  the  difference  between  the  experimental 
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Table  II. 

Averaged  Experimental/Theoretical  Amplitudes  for  the  Inlet  40"  from  Entrance  Plane  (SPL  in  dB) 


\  ^req. 

\(Hz) 

DegreeSv 

3,f  <t  \ 

300 

350 

400 

450 

500 

550 

600 

650 

700 

750 

0 

96.6 

93.8 

96.6 

96.2 

98.9 

99.1 

— 

104.2 

103.1 

108.1 

108.4 

BQ 

110.3 

110.8 

B 

109.0 

106.6 

B 

11.25 

95.9 

93.8 

96.1 

96.1 

99.5 

99.1 

103.5 

103.0 

108.6 

108.3 

B 

110.0 

110.7 

B 

B 

ES 

22.5 

95.3 

93.6 

96.3 

95.9 

99.1 

98.8 

102.4 

102.7 

109.1 

108.0 

111.0 

112.5 

110.6 

110.2 

IRQ 

107.3 

105.9 

B 

33.75 

94.7 

93.3 

97.0 

95.6 

98.1 

98.5 

B 

107.8 

107.4 

110.7 

111.9 

m 

107.4 

106.7 

B 

B 

45 

93.9 

93.0 

96.9 

95.2 

97.9 

98.0 

103.4 

101.7 

105.8 

106.7 

111.4 

111.1 

108.5 

108.6 

B 

104.6 

104.1 

B 

56.25 

92.7 

92.5 

95.8 

94.7 

1 

101.0 

101.1 

104.7 

106.0 

110.1 

110.1 

107.0 

107.5 

103.8 

104.5 

B 

B 

67.5 

91.1 

92.0 

95.0 

94.2 

98.7 

96.9 

100.1 

100.4 

m 

109.0 

109.1 

m 

103.0 

103.1 

mu 

78.75 

91.6 

91.3 

93.7 

93.5 

97.2 

96.2 

101.1 

99.7 

1 

107.0 

108.3 

m 

102.0 

101.9 

100.3 

99.8 

101.9 

98.8 

90 

91.0 

90.7 

92.3 

92.8 

95.9 

95.5 

101.5 

99.0 

a 

107.2 

107.6 

m 

100.7 

100.9 

B 

99.9 

97.4 

Driver 

127.0 

127.0 

127.0 

127.0 

127.0 

127.0 

127.0 

127.0 

127.0 

127.0 

'Vverage 

\bsolute 

irror 

1.18 

0.72 

0.54 

1.01 

0.54 

0.60 

0.59 

B 

1.09 

2.03 

Table  III 


Averaged  Experimental/Theoretical  Phases  for  the  Inlet  W  from 
the  Entrance  Plane  (Degrees  Relative  to  the  Driver) 


and  theoretical  results  in  prospective.  First,  there  are  the  obvious 
experimental  errors  caused  by  microphone  amplifier  drift,  temperature 
changes  in  the  anechoic  chamber  of  and  the  microphone  placement  in  the 
chamber.  Comparing  test  results  and  calibrations  these  errors  are  estimated 
to  account  for  up  to  0.5  dB  in  amplitude  and  5  degrees  in  phase.  Another 
source  of  error  is  the  anechoic  chamber  which,  as  stated  before,  can  account 
for  errors  up  to  3  dB  and  10  degrees.  A  more  subtle  source  of  experimental 
error  is  due  to  scale  switching  on  the  microphone  amplifiers  which  can 
account  for  as  much  as  0.5  dB  (i.e.  100  dB  on  the  90  to  110  dB  scale  reads  as 
100.5  dB  on  the  100  to  120  dB  scale). 

The  theoretical  errors  are  estimated  by  comparison  with  exact 
solutions  for  similar  geometries.  It  is  found  that  the  computer  programs 
introduce  about  1%  error,  insignificant  in  dB  but  accounting  for  as  much  as  5 
degrees  in  phase.  Another  source  of  error  is  the  assumption  of  a  plane  wave 
at  the  driver  plane;  unfortunately  the  effect  of  this  assumption  can  not  be 
easily  estimated.  As  none  of  the  above  mentioned  errors  are  specifically 
geometry  dependent,  these  error  estimates  apply  to  the  straight  duct 
configurations  too. 

The  results  of  the  experimental  tests  and  of  the  computer  runs  for 
the  hard  walled  straight  duct  configurations  are  presented  in  Tables  IV  and 
V.  The  amplitude  results  are  tabulated  in  decibels  in  Table  IV  while  the 
phase  results  in  degrees  appear  in  Table  V.  Again,  the  experimental  results 
presented  are  actually  the  average  of  two  separate  runs. 

It  can  be  seen  that  all  of  the  amplitude  results  are  very  good  as  the 
experimental  and  theoretical  values  at  each  test  frequency  have  average 


29 


Table  IV 


Averaged  Experimental/Theoretical  Amplitudes  for  the  Hard  Walled  Straight  Duct  40" 

from  the  Entrance  Plane  (SPL  in  dB) 


\  Freq . 

X.  (Hz) 
Degrees 

300 

350 

400 

450 

500 

550 

600 

650 

700 

0 

102.3 

99.8 

104.6 

104.2 

110.8 

110.5 

— 

114.1 

113.9 

110.1 

110.4 

B 

B 

B 

108.6 

108.1 

11.25 

1 

104.9 

104.1 

111.4 

110.4 

n 

110.9 

110.3 

107.0 

108.0 

106.4 

107.1 

B 

107.7 

107.8 

22.5 

a 

104.8 

103.8 

110.3 

110.0 

113.1 

113.4 

110.4 

109.8 

B 

B 

106.8 

106.5 

106.1 

107.1 

33.75 

n 

m 

110.0 

109.5 

113.1 

112.7 

108.1 

109.1 

u 

B 

B 

m 

45 

99.1 

98.8 

104.4 

102.9 

109.7 

108.8 

■ 

B 

107.2 

105.6 

I 

| 

m 

56.25 

97.9 

98.3 

H 

B 

■ 

106.3 

106.9 

B 

B 

m 

67.5 

96.9 

97.7 

102.3 

101.6 

108.3 

107.2 

B 

B 

H 

102.6 

101.5 

101.4 

101.0 

99.3 

101.1 

78.75 

96.7 

97.0 

101.5 

100.9 

107.0 

106.4 

110.2 

108.9 

B 

n 

101.5 

99.8 

99.9 

99.3 

100.5 

99.3 

90 

96.7 

96.3 

99.7 

100.2 

105.2 

105.8 

108.4 

108.1 

B 

100.7 

100.2 

98.9 

98.4 

98.8 

97.6 

95.5 

97.6 

Driver 

126.9 

126.9 

126.9 

126.9 

126.9 

126.9 

126.9 

126.9 

126.9 

Average 

Absolute 

irror 

0.96 

0.96 

0.67 

0.70 

0.71 

0.77 

0.91 

0.68 

1.00 
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absolute  errors  of  one  decibel  or  less.  The  phase  results  also  compare  very 
well  except  at  550  Hz  where  the  average  absolute  error  is  about  twice  that 
at  the  other  frequencies.  The  probable  reason  for  this  is  that  this  frequency 
is  very  close  to  the  tuning  frequency  of  the  liner,  which  is  about  558  Hz  (See 
Fig.  6.),  and  that  the  liner  orifices  were  closed  off  with  tape  for  the  hard 
walled  tests  which  may  not  have  been  100  percent  effective  in  keeping  the 
liner  from  influencing  the  phase  results.  It  should  be  noted,  however,  that 
this  method  of  closing  off  the  liner  was  effective  as  far  as  the  amplitude 
results  are  concerned  (See  Table  IV.).  The  amplitude  and  phase  results  at  550 
Hz  are  plotted  in  Fig.  9  of  Appendix  C  for  this  case. 

The  experimental  and  test  results  for  the  case  of  the  soft,  lined 
walled  straight  duct  configuration  are  tabulated  in  Tables  VI  and  VII  for  the 
amplitude  and  phase,  respectively.  The  results  for  the  amplitude  and  phase 
at  500  Hz  are  plotted  in  Fig.  10  of  Appendix  C.  The  errors  in  the  amplitude 
are  in  general  less  than  3  dB  and  the  errors  in  the  phase  less  than  10  degrees 
except  at  550  Hz  where  the  errors  for  both  the  amplitude  and  the  phase  are 
very  large.  These  differences  between  the  theoretical  and  experimental 
results  are  probably  due  to  the  liner  theory  used  for  calculating  the 
admittance  of  the  liner  and  machining  errors  in  the  manufacture  of  the 
liner. 
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Table  VI 

Averaged  Experimental/Theoretical  Amplitudes  for  the  Soft  Walled  Straight  Duct  40" 

from  the  Entrance  Plane  (SPL  in  dB) 


350  400  450  500  550  6 00  650  700 


Table  VII 

Averaged  Experimental/Theoretical  Phases  for  the  Soft  Walled  Straight  Duct  40" 
from  the  Entrance  Plane  (Degrees  Relative  to  the  Driver) 


\  Freq. 

(Hz) 

DegreeSv 

Oif  <t  \ 

300 

350 

400 

450 

500 
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600 

650 

700 

1 
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11 
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n 
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68 

■59 
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4 
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80 
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69 

KB 
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9 
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-45 

16 
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76 
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22.5 
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4 

134 
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-59 

82 
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68 
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-71 

12 
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-47 

9 

-57 

75 

Bl 

33.75 

-61 

5 
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-128 
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89 
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68 

■a 
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n 
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m 

-55 

n 
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62 
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45 
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H 
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Ba 
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El 
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67 
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mm. 
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-52 

■a 
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72 

156 

56.25 

wm 
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B9 
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67 
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-61 

mm 

m 

-43 

6 
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85 
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67.5 

-56 

U 

wm 

Ba 

-54 
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69 
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wm 

wm 

-120 

-36 

19 

-48 

69 
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78.75 

m 

m 

-116 

-49 
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72 

145 
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25 

wm 
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bi 
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86 

B9 

90 
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19 

mm 
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H 
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80 

Ea 
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■ 

0 

0 

0 

0 

0 

0 

0 

0 
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■■ 

m 

Absolute 

4.8 

9.3 

3.8 
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8.1 

6.1 

illll 

Error 

m 

HI 

■H 
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F.  Determination  of  the  "Effective"  Liner  Admittance 


As  stated  in  the  previous  sections  normal  machining  errors  can 
change  the  resonant  frequency  of  the  liner  by  as  much  as  1.5  Hz  which  can 
change  the  sound  field  by  as  much  as  3  dB  and  5  degrees  due  to  the  liners 
"peaky"  absorption  curve  (See  Fig.  6.).  Even  this  can't,  however,  account  for 
the  large  differences  between  the  experimental  and  theoretical  results  at 
550  Hz  for  the  lined  straight  duct.  Noteing  that  the  calculated  results  were 
consistently  lower  than  those  measured,  it  was  assumed  that  the  liner 
theory^  predicted  the  admittance  (i.e.  the  effectiveness)  of  the  liner  to  be 
higher  than  it  actually  was.  To  check  this  hypothesis  a  parametric  set  of 
computer  runs  was  conducted  in  which  the  admittance  was  systematically 
varied  to  see  if  the  true  "effective"  admittance  of  the  liner  at  550  Hz  could 
be  determined  by  comparison  with  the  experimental  results.  From  this  set  of 
computer  runs  it  was  found  that  the  true  "effective"  admittance  of  the  liner 
was  approximately 


y  =- 0.13- i0.28  @  550  Hz 

(10) 

which  is  the  real  part  of  the  admittance  calculated  from  Eqs.  (2-4)  and  (8) 
multiplied  by  0.35  and  the  imaginary  part  multiplied  by  0.65.  The 
experimental  and  theoretical  results,  calculated  using  the  admittance  in  Eq. 
(10),  are  tabulated  in  Table  VIII.  As  can  be  seen  these  results  are  more  in 
line  with  what  was  found  at  the  other  frequencies  (See  Tables  VI  and  VII.),  so 
it  is  concluded  that  the  Garrison  theory  of  Ref.  (5)  does  not  give  the  proper 
admittance  for  this  type  of  liner  near  resonance  (i.e.  it  predicts  the 


4 


Table  VIII 

Averaged  Experimental/Theoretical  Acoustic  Radiation  40"  from  the 
Entrance  Plane  of  the  Soft  Walled  Straight  Duct  at  550  Hz  with 

y  =-  0.13  -  i  0.28 


admittance  to  be  too  high).  This  exercise  also  points  out  that  the  integral 
equation  technique  can  be  used  to  determine  a  liners  true  "effective" 
admittance  from  experimental  measurements  in  the  field. 


37 


r\n.: 


« 


References 

1. )  Meyer,  W.  L.  ,  Bell,  W.  A.,  and  Zinn,  B.  T.,  "Prediction  of  the  Sound 

Field  Radiated  from  Axisymmetric  Surfaces,"  AIAA  Paper  No.  78- 
195,  presented  at  the  AIAA  16th  Aerospace  Sciences  Meeting, 
Huntsville,  Alabama,  January  16-18,  1978. 

2. )  Meyer,  W.  L.,  Bell,  W.  A.,  Stallybrass,  M.  P.,  and  Zinn,  B.  T., 

"Prediction  of  the  Sound  Field  Radiated  from  Axisymmetric 
Surfaces,"  Journal  of  the  Acoustical  Society  of  America,  Vol.  63,  No. 
2,  pp.  631-638,  March  1979. 

3. )  Meyer,  W.  L.,  Bell,  W.  A.,  and  Zinn,  B.  T.,  "Sound  Radiation  from 

Finite  Length  Axisymmetric  Ducts  and  Engine  Inlets,"  AIAA  Paper 
No.  79-0675,  presented  at  the  AIAA  5th  Aeroacoustics  Conference, 
Seattle,  Washington,  March  12-14,  1979. 

4. )  Miller,  B.  A.,  Dastoli,  B.  J.,  and  Wesoky,  H.  L.,  "Effect  of  Entry-Lip 

Design  on  Aerodynamics  and  Acoustics  of  High-Throat-Mach-Number 
Inlets  for  the  Quiet,  Clean,  Short-Haul  Experimental  Engine,"  NASA 
TM  X-3222,  Lewis  Research  Center,  Cleveland,  Ohio,  May  1975. 

5. )  Garrison,  G.  D.,  "A  Study  ofthe  Suppression  of  Combustion 

Oscillations  with  Mechanical  Damping  Devices,"  Phase  II  Summary 
Report,  Pratt  &  Whitney  Aircraft  Co.,  PWA  FR-1922,  July  15,  1966. 


« 


38 


Appendix  A 

Publications  Generated  by  this  Research  Effort 

A. )  "Integral  Solution  of  Three  Dimensional  Acoustic  Radiation 

Problems,"  W.  L.  Meyer,  W.  A.  Bell  and  B.  T.  Zinn,  Proceedings  of  the 
International  Symposium  on  Innovative  Numerical  Analysis  in  Applied 
Engineering  Science,  Versailles,  France,  May  23-27,  1977. 

B. )  "Predicting  the  Acoustics  of  Arbitrarily  Shaped  Bodies  Using  an 

Integral  Approach,"  W.  A.  Bell,  W.  L.  Meyer  and  B.  T.  Zinn,  AIAA 
Journal,  Vol.  15,  No.  6,  pp.  813-820,  June  1977. 

C. )  "Prediction  of  the  Sound  Field  Radiated  from  Axisymmetric 

Surfaces,"  W.  L.  Meyer,  W.  A.  Bell  and  B.  T.  Zinn,  AIAA  Paper  No. 
78-195,  presented  at  the  AIAA  16th  Aerospace  Sciences  Meeting, 
Huntsville,  Alabama,  January  16-18,  1978. 

D. )  "Boundary  Integral  Solutions  of  Three  Dimensional  Acoustic  Radiation 

Problems,"  W.  L.  Meyer,  W.  A.  Bell,  M.  P.  Stallybrass  and  B.  T.  Zinn, 
Journal  of  Sound  and  Vibration,  Vol.  59,  No.  2,  pp.  245-262,  July  1978. 

E. )  "Prediction  of  the  Sound  Field  Radiated  from  Axisymmetric 

Surfaces,"  W.  L.  Meyer,  W.  A.  Bell,  M.  P.  Stallybrass  and  B.  T.  Zinn, 
Journal  of  the  Acoustical  Society  of  America,  Vol.  63,  No.  2,  pp.  631- 
638,  March  1979. 


39 


F. )  "Sound  Radiation  from  Finite  Length  Axisymmetric  Ducts  and  Engine 

Inlets,"  W.  L.  Meyer,  W.  A.  Bell  and  B.  T.  Zinn,  AIAA  Paper  No.  79- 
0675,  presented  at  the  AIAA  5th  Aeroacoustics  Conference,  Seattle, 
Washington,  March  12-14,  1979. 

G. )  "Acoustic  Radiation  from  Axisymmetric  Ducts:  A  Comparison  of 

Theory  and  Experiment,"  W.  L.  Meyer,  B.  R.  Daniel  and  B.  T.  Zinn, 
AIAA  Paper  No.  80-0097.  presented  at  the  AIAA  18th  Aerospace 
Sciences  Meeting,  Pasadena,  California,  January  14-16,  1980. 


Appendix  B 

Presentations  Generated  by  this  Research  Effort 
"Integral  Solutions  of  Three  Dimensional  Acoustic  Radiation 
Problems,"  Presented  at  the  International  Symposium  on  Innovative 
Numerical  Analysis  in  Applied  Engineering  Science  at  Versailles, 
France,  May  23-27,  1977.  (Zinn) 

"Prediction  of  the  Sound  Field  Radiated  from  Axisymmetric 
Surfaces,"  Presented  at  the  AIAA  16th  Aerospace  Sciences  Meeting 
at  Huntsville,  Alabama,  January  16-18,  1978.  (Bell) 

"Sound  Radiation  from  Finite  Length  Axisymmetric  Ducts  and  Engine 
Inlets,"  Presented  at  the  5th  AIAA  Aeroacoustics  Conference  at 
Seattle,  Washington,  March  12-14,  1979.  (Zinn) 

"Acoustic  Radiation  from  Axisymmetric  Ducts:  A  Comparison  of 
Theory  and  Experiment,"  Presented  at  the  AIAA  18th  Aerospace 
Sciences  Meeting,  Pasadena,  California,  January  14-16,  1980.  (Meyer) 


Appendix  C 


AIAA  Paper  Number  80-0097  Presented  at  the  AIAA  18th  Aerospace 
Sciences  Meeting,  Pasadena,  Calif.,  3an.  14-16,  1980. 


42 


AIAA-80-0097 


Acoustic  Radiation  from 


Axisymmetric  Ducts:  A  Comparison 
of  Theory  and  Experiment 


W.  L.  Meyer,  B.  R.  Daniel  and 
B.  T.  Zinn,  Georgia  Institute  of 
Technology,  Atlanta,  Ga. 


Al  AA  1 8th 

AEROSPACE  SCIENCES  MEETING 


January  14-16,  1980/Pasadena,  California 


fm  parmlaaiaa  to copy  or  rapaWtih,  contort  tho  Amortcan  Imtltuto  ot  Aoronouttci  and  Astronautici>1290  Avonut  at  tha  Amarieai,  Now  York,  N  Y.  10019 


ACOUSTIC  RADIATION  FROM  AXISYMMETRIC  DUCTS: 
A  COMPARISON  OF  THEORY  AND  EXPERIMENT* 

W.  L.  Meyer  ,  B.  R.  Daniel*  and  B.  T.  Zinn** 

School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Ga.  30332 


Abstract 

A  special  integral  representation  of  the  exterior 
solutions  of  the  Helmholtz  equation  is  used  to  calculate  the 
free  field  acoustic  radiation  patterns  around  two  finite 
axisymmetric  bodies;  a  straight  pipe  and  a  jet  engine  inlet. 
The  radiation  patterns  around  these  bodies  are  then 
measured  experimentally,  with  the  free  field  being 
approximated  through  the  use  of  an  anechoic  chamber.  The 
inlet  tested  has  a  hard  wall  while  the  straight  pipe  is  tested 
with  both  a  hard  and  a  lined  wall.  The  computed  theoretical 
and  the  measured  experimental  acoustic  radiation  patterns 
are  found  to  be  in  good  agreement.  A  discussion  of  possible 
sources  of  error,  both  theoretical  and  experimental,  is 
included. 

Introduction 

One  of  the  major  problems  facing  the  aircraft 
industry  today  is  to  reduce  the  noise  radiated  to  the  ground 
from  aircraft  engines  without  sacrificing  any  of  the  overall 
efficiency  of  the  aircraft.  A  major  source  of  engine  noise  is 
the  compressor  or  the  fan  noise  which  is  radiated  out  of  the 
jet  engine  through  its  inlet  section.  Research  efforts 
directed  toward  reducing  these  noise  sources  have  included, 
among  other  things,  the  reduction  of  the  sound  level  in  the 
jet  engine  inlet  section,  by  adding  sound  absorbing 
materials  (i.e.,  acoustic  liners)  to  the  inside  of  the  engine 
inlet.  Weight  and  volume  of  these  acoustic  liners  are  of 
prime  concern  to  the  aircraft  industry  as  these  are  directly 
related  to  the  overall  efficiency  of  the  aircraft;  thus, 
efficient  acoustic  liner  designs  are  sought.  In  the  past, 
optimum  liner  designs  have  been  found  by  extensive  full 
scale  testing  of  engines  with  various  liner  configurations 
which  is  a  very  costly  and  time  consuming  process.  Most  of 
this  testing  could  be  eliminated  by  the  development  and  use 
of  efficient,  accurate  analytical  procedures  for  the 
prediction  of  the  sound  field  radiated  from  lined  jet  engine 
inlets.  The  development  of  such  a  method,  based  upon  a 
special  integi  al  representation  of  the  radiation  solutions  of 
the  Helmholtz  equation,  has  been  discussed  at  length  by  the 
authors  of  this  paper  in  earlier  publications'1'  .  In  the 
present  study  the  applicability  of  this  integral  solution 
technique  is  investigated  by  comparing  its  predictions  with 
the  results  obtained  from  an  experimental  study. 

Background 

In  this  paper  a  special  axisymmetric  integral 
representation  of  the  exterior  solutions  of  the  Helmholtz 
equation  u’4'  is  employed  to  theoretically  calculate  the 
free  field  acoustic  radiation  patterns  surrounding  two  finite 
axisymmetric  bodies;  namely,  a  straight  duct  and  a  jet 
engine  inlet  configuration.  The  inlet  used  in  these  studies  is 
the  so  called  QCSEE  inlet  (i.e.,  the  quiet,  clean,  short-haul 
experimental  engine  inlet  of  Reference  5).  In  previous 
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studies  by  this  group  the  sound  fields  radiated  from 
complicated  geometries  with  complex  (mixed)  boundary 
conditions  were  calculated  using  a  solution  procedure  based 
on  the  aforementioned  integral  representation.  The 
solutions  generated  were  found  to  be  in  excellent 
agreement  with  "exact"  solutions  calculated  by  employing 
the  method  of  Separation  of  Variables.  A  detailed 
development  of  the  integral  equations  and  solution 
procedures  used  along  with  many  comparisons  with  exact 
solutions  for  various  bodies,  both  3-D  and  axisymmetric, 
can  be  found  in  Refs.  (1-4  and  6).  In  this  connection  it 
should  be  pointed  out  that  this  theoretical  technique  and 
solution  procedure  have  been  found  to  be  both  accurate  and 
computationally  efficient  when  compared  to  other 
methods. 

Theoretical  Method 

In  this  study  a  so  called  "integral  technique"  is 
employed  to  calculate  the  sound  radiated  from  various 
axisymmetric  configurations.  The  particular  method  is 
unique  in  that  it  is  applicable  at  all  nondimensional  wave 
numbers  ka  (where  k  is  the  wave  number  and  a  is  ar 
appropriate  body  dimension),  it  contains  no  tangential 
derivatives  on  the  surface  of  the  body  S  (See  Fig.  1.)  and  it 
contains  no  singular  kernels  which  cannot  be  handled 
numerically  by  straight  forward  means. 


Figure  1.  Definition  of  Elements  on  the  Body. 


The  details  of  the  derivation  of  the  particular  integral 
equations  used  are  presented  in  Refs.  (1,  2  and  6)  and 
therefore  will  not  be  repeated  here.  The  basic  integral 
equation  employed  on  the  surface  of  the  body  is  written  in 
terms  of  the  acoustic  potential  tp  and  it  contains  surface 
integrals  over  the  surface  of  the  body  S.  If  we  now  define 
fras  the  outward  normal  derivative  from  the  surface,  P 
and  Q  as  points  on  the  surface  with  their  unit  outward 
normals  defined  as  and  ^  respectively,  and  as  the 
normal  acoustic  velocity  this  integral  equation  takes  the 
following  general  form 
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In  this  equation  G(P,Q)  can  be  any  fundamental  solution  of 
the  Helmholtz  equation.  In  this  study  it  has  been  chosen  to 
be  the  free  space  Green's  Function  of  the  Helmholtz 
equation;  that  is 
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Figure  2.  General  Axisymmetric  Geometry  and 
2-D  Projection. 


•vhere  r(P,Q)  is  the  distance  between  the  points  P  and  Q  on 
the  surface  of  the  body. 


Once  the  acoustic  potential  is  calculated  on  the 
surface  of  the  body  using  Eqn.  (1),  the  following  integral 
representation  for  tp  may  be  used  to  calculate  the 
acoustic  potential  anywhere  in  the  field  surrounding  the 
body.  In  this  case  the  point  P  is  no  longer  on  the  body  S  but 
in  the  field  surrounding  it. 
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If  the  body  of  interest  is  axisymmetric,  as  most  vat 
engine  inlets  are,  certain  simplifications  may  be  made'3’  . 
Let  us  now  define  a  new  cylindricaily  symmetric  acoustic 
potential  and  normal  acoustic  velocity  as 


and  s  is  the  distance  along  the  body  in  the  o-z  plane  (e.g. 
0  <  s  <t).  In  doing  this  we  have  effectively  separated  out 
all  the  tangential  modes;  thus,  the  formulation  that  will  be 
presented  for  axisymmetric  bodies  will  be  cylindricaily 
symmetric  in  that  all  the  tangential  modes  (i.e.,  m  - 
0,1,2...)  may  be  solved  for  separately.  Having  these 
definitions,  Eqn.  (1)  can  be  rewritten  as 
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where  s  is  the  distance  along  the  2-D  projection  of  the 
body  in  the  p-z  plane  (See  Fig.  2.).  The  kernel  functions 
Kj  and  K2  are  defined  as 
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where  m  is  the  tangential  mode  number  and  6  is  defined  as 
in  Figure  2 
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in  the  field  suriounding  the  Body,  inese  were  then 
compared  with  the  acoustic  potential  that  the  simple 
source  would  generate  at  each  point  in  the  field.  Very  good 
agreement  was  found;  for  53  integration  points  on  the  body 
along  s  (See  Fig.  2.)  the  error  was  always  below  10%  in 
the  real  and  imaginary  parts  of  4  for  nondimensional  wave 
numbers  up  to  ka  *  10.  It  was  also  found  that  increasing 
the  number  of  integration  points  generally  decreased  the 
error  proportionally  (e.g.  using  102  points  the  maximum 
error  in  the  acoustic  potential  at  a  field  point  was  about 
5%  for  the  inlet  at  a  nondimensional  wave  number  of  ka  = 
10.).  The  results  and  details  of  many  such  tests  are 
presented  in  Ref.(3). 


and  the  influence  functions  are  defined  as 
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The  Test  Bodies 

The  integral  formulation  has  been  used  in  the 
present  study  to  calculate  the  sound  fields  radiated  from 
two  axisymmetric  bodies,  a  straight  duct  and  a  jet  engine 
inlet.  The  theoretical  and  experimental  configurations  do 
not  compare  exactly  as  an  accurate  description  of  the 
"back  side"  of  the  experimental  bodies  would  entail  the  use 
of  too  many  theoretical  points  on  the  bodies.  Therefore, 
different  external  terminations  were  given  to  the 
theoretical  models  in  the  interest  of  conserving  computing 
time  and  computer  storage  space.  In  this  connection  it  has 
been  found  through  theoretical  studies,  that  the  exact  form 
of  the  rear  termination  of  the  body  has  little  effect  on  the 
sound  field  radiated  in  the  forward  half  plane.  Since  we  are 
only  really  interested  in  the  sound  field  in  the  forward  half 
plane,  this  approximation  is  not  considered  to  be  a  major 
source  of  error  in  this  study. 


One  should  note  that  8p  has  been  chosen  to  be  zero  (i.e.  The  first  axisymmetric  body  tested  was  a  straight 

cos  m0  =1)  in  the  above  equation  which  can  be  done  duct-  Beiow  is  a  sketch  of  the  theoretical  model  employed, 
without  Sny  loss  of  generality. 


Equation  (3)  may  also  be  rewritten  in  the  form 
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As  can  be  seen,  this  formulation  of  the  problem  reduces  the 
solutions  of  Eqns.  (1)  and  (3)  to  the  evaluation  of  line 
Integrals  on  the  2-D  projections  of  the  body  (See  Fig.  2.). 
Also,  this  formulation  can  account  for  tangential  modes; 
which  must,  however,  be  solved  for  separately. 

The  theoretical  method  described  above  was 
checked  for  accuracy  by  employing  several  axisymmetric 
geometries  including  a  straight  duct  and  a  jet  engine  inlet. 
The  calculated  radiation  patterns  generated  by  this  method 
were  compared  with  various  "exact"  solutions.  These  exact 
solutions  were  found  by  assuming  that  some  simple  sources 
(e.g.  monopoles,  dipoles,  and  quadrupoles)  were  located 
within  the  body  and  then  calculating  the  normal  acoustic 
velocity  V  and/or  the  admittance,  defined  as  y«v/»,at 
points  on  the  surface  of  the  body.  These  values  were  then 
used  as  the  boundary  conditions  on  the  surface  of  the  body 
and  the  acoustic  potential  was  calculated  at  various  points 


Figure  3.  Theoretical  Straight  Duct  Model. 


Its  L/a  equals  2.1 10,  the  same  as  the  experimental  model. 
It  will  be  noted  that  the  theoretical  model  has  a 
hemispherical  termination  which  differs  from  the  rear 
termination  of  the  experimental  model  shown  below. 
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Figure  4.  Experimental  Straight  Duct  Test 
Configuration. 


In  the  experimental  model  the  driver  is  placed  at  the  throat 
of  a  nozzle  section,  which  assures  that  the  sound  waves  are 
plane  at  the  driver  plane  where  a  1/4"  condensor 
microphone  is  located  to  provide  a  reference  pressure  level 
in  decibels.  All  of  the  tests  were  conducted  at 
nondimensional  wave  numbers  ka  below  the  l.T  mode  of  the 
duct  (i.e.  k  <  1.84)  and  thus  we  were  assured  of  a  plane 
wave  at  the*  driver  plane.  This  has  been  checked  out 
experimentally  by  sweeping  the  reference  microphone 
radially  across  the  duct  and  it  has  been  found  to  be  true 
within  1  dB  in  amplitude  and  5°  in  phase  (i.e.  there  is  less 
than  1  dB  and  5°  variation  between  the  wall  and  the  center 
of  the  duct).  In  the  theoretical  model  the  existence  of 
plane  wave  excitation  at  the  driver  plane  is  assumed. 

The  straight  duct  was  tested  in  two  configurations; 
that  is,  hard  wall  and  lined  wall  configurations.  The  lined 
wall  configuration  consisted  of  180  Helmholtz  resonators  (9 
axial  rows  by  20  radial  rows).  A  sketch  of  one  of  the 
resonators  is  shown  below. 


Figure  5.  Helmholtz  Resonator. 


For  the  hard  walled  tests  the  small  holes  inside  the  duct 
were  simply  covered  with  tape. 

The  second  axisymmfttric  body  tested  was  a  model 
of  an  actual  engine  inlet'  .  It  has  an  L/a  of  2.0.  The 
theoretical  model  has  a  hemispherical  termination  similar 
to  the  one  used  for  the  straight  duct.  A  sketch  of  the 
theoretical  inlet  model  is  presented  below. 


Figure  6.  Theoretical  Inlet  Model. 


In  the  experimental  model  the  centerbody  was  held  in  place 
by  4  small  wing  cross-sectional  struts  set  at  90°  angles. 
The  microphone  measurement  for  the  amplitude  and  phase 
at  the  driver  plane  was  made  half  way  between  2  of  these 
struts.  Again,  the  small  condensor  microphone  was  swept 
across  the  driver  plane  radially  to  check  for  the  presence 
of  plane  wave  excitation  which  was  found  to  exist  within 
the  same  limits  as  for  the  straight  duct  (i.e.,  less  than  a  1 
dB  change  in  amplitude  and  5°  in  phase  across  the  driver 
plane).  A  sketch  of  the  experimental  inlet  model  is 
presented  below. 
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Figure  7.  Experimental  Inlet  Test  Configuration. 


» 


Again,  the  exact  mathematical  form  of  the  curves  that 
make  up  the  inlet  contours  can  be  found  in  Reference  5. 

Results 

In  the  anachoic  chamber  tests  conducted  under  this 
program,  the  field  measurements  were  taken  40"  from  the 
center  of  the  duct  exit  plane  on  a  circular  arc  at  1 1  1/4° 
increments  from  the  centerline  of  the  duct.  One-half  inch 
5  &  K  microphones  were  used  in  these  tests  for  the  field 
measurements.  The  anechoic  chamber  used  in  the 
experiments  has  interior  dimensions  of  10'  x  13'  x  6  1/2* 
high.  The  acoustic  insulation  used  in  the  chamber  is 
fiberglass  which  is  approximately  2'  thick.  The  ducts  tested 
(i.e.,  the  straight  duct  and  the  inlet)  and  the  microphones 
used  for  the  field  measurements  were  3*  off  the  floor.  A 
plan  view  of  the  anechoic  chamber  with  a  typical  test 
set-up  is  presented  below. 


Figure  8.  Anechoic  Chamber  (Plan  View  )  with 
Test  Set-Up. 


Free  field  measurements  were  made  in  the  chamber 
to  check  its  "anechoicness".  For  the  distances  of  interest  in 
the  chamber  (i.e.  up  to  4(7'  from  a  source)  free  field 
conditions  were  generally  approximated  to  within  3  dB  in 
amplitude  and  10°  in  phase. 

The  procedure  for  the  tests  was  to  take  the 
measured  SPL  in  dB  from  the  1/4"  B  it  K  microphone  at  the 
driver  plane  and  input  the  data  into  the  theoretical  model 
assuming  plane  wave  excitation  at  the  driver  plane.  Again, 
this  was  checked  experimentally  for  each  configuration  and 
found  to  be  generally  true  within  1  dB  in  amplitude  and  5 
degrees  in  phase.  Since  a  plane  wave  was  assumed,  the 
phase  was  not  measured  at  the  driver  plane  and  zero  phase 
was  assumed  for  this  location  in  the  theoretical  model.  This 
can  be  done  without  any  loss  of  generality  as  the  phase 
differences  between  the  field  microphones  and  the  driver 
plane  microphone  are  the  quantities  of  interest  which  were 
measured  in  this  program.  Once  the  driver  plane  acoustic 
pressure  was  input  into  the  theoretical  model,  the  far  field 
sound  distribution  was  calculated  and  compared  with  the 
measured  experimental  data. 

The  computer  used  for  these  analyses  was  the 
Georgia  Tech  Cyber  70/74.  Typical  run  times  to  calculate 
the  distribution  of  the  acoustic  potential  c  (i.e..  See  Eqn. 
(5).)  with  100  points  on  the  theoretical  body  are  about  6 
minutes  and  to  calculate  the  acoustic  potential  in  the  far 
field  at  20  points  (i.e.,  See  Eqn.  (9).)  are  about  1  minute.  In 
the  theoretical  calculations  of  the  surface  potential  on  the 
straight  duct  102  points  were  used  while  in  the  calculations 
for  the  inlet  (See  Fig.  5.)  97  points  were  used. 

In  the  numerical  integration  oi  Eqns.  (5)  and  (9)  a  2 
point  integration  formula  was  used  in  the  s  direction  (See 
Fig.  2.).  It  should  be  noted,  however,  that  an  even 
integration  formula  must  be  used  here  as  an  odd  formula 
would  place  a  point  in  the  center  of  the  integration  region 
which  would  cause  the  various  kernel  functions  (i.e.,  See 
Eqns.  (6)-(8).)  to  go  to  infinity  when  r(P,Q)  goes  to  zero  as 
the  points  P  and  Q  coincide.  In  all  the  investigated  cases  a 
96  point  Gaussian  integration  formula  was  used  in  the 
0  direction. 

Tests  with  the  straight  duct  configuration  with  both  hard 
and  soft  walls  were  conducted  in  the  frequency  range  of 
300  to  700  Hz  at  50  Hz  increments.  Over  this  frequency 
range  a  reasonably  plane  wave  could  be  excited  at  the 
driver  plane*  .  Comparisons  of  calculated  and  measured 
data  for  the  hard  walled  straight  duct  configuration  are 


In  this  study  the  lower  frequency  limit  was  imposed 
by  the  limitations  of  the  University  driver. 
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presented  in  Fig.  9  for  550  Hz.  The  SPL  results  at  this 
frequency  are  representative  of  those  at  all  other 
frequencies  while  the  phase  results  are  the  worst  at  this 
frequency.  The  probable  reason  for  this  is  that  550  Hz  is 
very  close  to  the  tuning  frequency  of  the  liner,  which  was 
calculated  to  be  around  558  Hz  under  these  test  conditions 
and  the  tape  used  to  close  off  the  liner  holes  for  the  hard 
walled  tests  was  not  100  percent  effective.  The  average 
absolute  errors  for  the  amplitude  in  dB  and  for  the  phase  in 
degrees  are  presented  in  Table  I  for  all  of  the  tests  run. 


550  HZ 


POSITION 

(DEGREES  OFF  DUCT 


Figure  9.  Acoustic  Radiation  40"  from  the 
Hard  Walled  Straight  Duct. 


Comparisons  of  experimental  and  calculated  far 

field  pressures  for  the  lined,  straight  wall  duct  are 
presented  in  Fig.  10  at  a  frequency  of  500  Hz.  These  results 
are  typical  of  those  obtained  at  all  but  one  of  the  other 
frequencies  tested.  At  550  Hz,  which  is  very  close  to  the 
calculated  resonant  frequency  of  the  liner,  the  results  show 
significantly  more  error  (See  Table  I.).  Since  the 

calculated  results  were  lower  than  those  measured  it  was 
assumed  that  the  liner  theory  results  predicted  the 

effectiveness  of  the  liner  to  be  higher  than  it  actually  was 
in  practice.  To  see  if  this  was  actually  the  case  a  set  of 
systematic  computer  runs  were  made  in  which  the 

effectiveness  of  the  liner  (i.e.,  its  admittance)  was 
reduced. 


CALCULATED  -  MEASURED  O 


Figure  10.  Acoustic  Radiation  40"  from  the 
Soft  WalJed  Straight  Duct. 


To  predict  the  sound  field  radiated  by  the  lined, 
straight  wall  duct  the  admittance, at  the  wall  is  required. 
In  this  case  available  liner  theory'"  was  used  to  predict  the 
admittance  of  the  Helmholtz  resonator  array.  For  the 
particular  resonators  used  (See  Fig.  5.)  in  this  study  the 
resonant  frequency  of  the  array,  fQ  ,  was  calculated  to  be 
around  558  Hz  and  the  specific  impedance  of  the  liner 
Z  =  Q  -  A  Z  was  found  to  be  given  by: 


KJ  -■  O.  0/39  8  V/ 

X  -  /o.  /  ae  (  ~j  ) 
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It  was  found  that  if  the  predicted  real  part  of  the 
admittance  was  taken  to  be  35  percent  effective  and  the 
imaginary  part  of  the  admittance  was  taken  to  be  65 
percent  effective  (i.e.  y  =  -  0.13  -  i  0.28  at  550  Hz)  than 
the  results  from  the  theoretical  calculations  were 
significantly  closer  to  the  experimentally  measured  data. 
The  average  absolute  errors  for  the  amplitude  in  dB  and  for 
the  phase  in  degrees  are  now  2.51  dB  and  5.00  degrees 
respectively  as  compared  to  the  much  larger  errors' 
incurred  using  the  admittance  values  predicted  by  the  liner 
theory  (See  Table  I.). 

Tests  for  the  inlet  configuration  were  also 
conducted  over  the  frequency  range  300  -  700  Hz  at  50  Hz 
increments.  Since  the  reference  lengths,  a,  are  slightly 
different  (See  Figs.  3,4,6  and  7.)  the  nondimensional  wave 
numbers  ka  are  different.  Comparisons  of  theoretical  and 
experimental  data  are  presented  in  Fig.  11  for  the  inlet 
configuration  at  600  Hz.  These  results  are  representative 
of  those  at  other  frequencies. 
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Discussion  of  Results 


600  Hz 
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Figure  11.  Acoustic  Radiation  90'  from  the  Inlet. 


The  causes  for  the  errors  appearing  in  the 
comparisons  of  the  last  section  will  be  briefly  discussed 
herein.  First,  there  are  the  obvious  experimental  errors 
caused  by  microphone  amplifier  drift,  temperature  changes 
in  the  anechoic  chamber  and  microphone  placement  in  the 
chamber.  These  errors  can  be  estimated  as  identical  tests 
were  run  on  different  days,  the  microphones  were 
calibrated  three  times  during  the  course  of  each  test,  and 
the  microphones  were  moved  and  reset  during  each  test  as 
there  are  9  positions  in  the  field  and  only  5  microphones 
were  used  in  each  test.  Comparing  test  results  and 
calibrations,  these  errors  are  estimated  as  being  of  the 
order  of  0.5  dB  in  amplitude  and  5  degrees  in  phase. 
Another  source  of  experimental  error  is  the  anechoic 
chamber  itself  which,  as  stated  before,  can  account  for  up 
to  3  dB  errors  in  amplitude  and  10°  errors  in  phase.  A 
more  subtle  source  of  error  in  the  lined  duct  tests  is  the 
liner  itself.  It  was  found  that  by  changing  the  dimensions  of 
the  Helmholtz  resonators  by  as  little  as  0.0005  inches  (i.e. 
common  machining  errors)  the  resonant  frequency 
calculated  for  the  liner  changed  by  as  much  as  1.5  Hz 
which  is  significant  for  this  type  of  liner  near  resonance 
due  to  its  highly  peaked  absorption  curve.  A  sample 
calculation  was  run  at  550  Hz  with  this  change  and  the 
calculated  results  changed  by  about  3  dB  in  amplitude  and  5 
degrees  in  phase.  Another  source  of  possible  error  for  this 
particular  case  is  obviously  the  imperfection  of  the  liner 
theory  itself;  the  determination  of  the  errors  caused  by  its 
shortcomings  are,  however,  beyond  our  current  capabilities. 

The  computer  programs  also  introduce  some  errors 
which  are  estimated  to  be  about  1  percent  by  comparing 
these  computer  results  with  exact  solutions  for  similar 
geometries  and  wave  numbers.  Although  these  errors  are 
insignificant  when  evaluated  in  dB,  they  can  be  as  high  as  5 
degrees  in  phase.  Another  source  of  error  is  the  assumption 
of  a  plane  wave  at  the  driver  plane;  the  effect  of  this  error 
or  the  results  in  the  far  field  cannot,  however,  be  easily 
estimated.  Other  sources  of  error  include  the  differences 
between  the  experimental  and  theoretical  geometries 
which  include  not  only  the  different  terminations  on  the 
back  side  of  the  bodies  but  also  the  stand  required  to  hold 
up  the  experimental  set-up  in  the  anechoic  chamber.  The 
errors  caused  by  these  differences  are  hopefully  small. 


300 

350 

900 

950 

500 

550 

600 

650 

700 

Hard 

Vailed 

Straight 

Duct 

dB 

0.96 

0.95 

0.66 

0.70 

0.71 

0.78 

0.92 

0.67 

1.02 

Deg. 

3.83 

5.56 

5.89 

9.61 

7.99 

13.78 

9.33 

7.72 

7.83 

Soft 

Vailed 

Straight 

Duct 

dB 

1.06 

1.96 

0.82 

0.63 

0.79 

15.71 

1.61 

0.92 

1.19 

Deg. 

9.67 

9.78 

9.28 

3.67 

9.33 

88.89 

8.11 

5.89 

6.61 

Inlet 

dB 

1.18 

0.72 

0.61 

1.08 

0.69 

0.63 

0.61 

0.38 

1.07 

Deg. 

3.17 

6.39 

3.83 

10.22 

5.00 

6.22 

9.00 

9.50 

5.00 

Table  I.  Average  Absolute  Errors. 
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Conclusions 


Acoustic  measurements  were  made  of  the  sound 
field  radiated  from  a  straight  duct  with  both  acoustically 
hard  and  soft  walls  and  a  jet  engine  inlet.  These 
measurements  were  then  compared  with  the  results  of  an 
integral  representation  of  the  solutions  of  the  Helmholtz 
equation  and  good  agreement  between  the  theoretical  and 
experimental  results  was  observed.  This  indicates  that  the 
integral  equations  used  and  the  techniques  employed  for 
solving  them  are  good  approximations  to  the  actual 
acoustic  behavior  of  arbitrarily  shaped  axisymmetric  ducts 
radiating  into  a  free  space.  This  is  significant  in  that  most 
theories  can  not  adequately  model  the  coupling  between 
the  acoustic  fields  inside  and  outside  a  duct.  Thus,  this 
technique  can  be  used  with  confidence  to  efficiently 
predict  the  sound  field  radiated  from  complex 
axisymmetric  geometries. 
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